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Abstract 



The constraint of a progressive decrease in residual renormalization scale dependence with increasing loop 
. order is developed as a method for obtaining bounds on unknown higher-order perturbative corrections to 

^ ' renormalization-group invariant quantities. This technique is applied to the inclusive semileptonic process 

, b uh'il' (explicitly known to two-loop order) in order to obtain bounds on the three- and four-loop pertur- 

■ bative contributions that are not accessible via the renormalization group. Combining this technique with the 

04 ' principle of minimal sensitivity, we obtain an estimate for the perturbative contributions to F (6 ^ uvel~) that 

incorporates theoretical uncertainty from as-yet-undetermined higher order QCD corrections. 

O 

' A variety of techniques have been developed for obtaining estimates of higher- loop corrections in perturbative QCD, 
, including the principle of minimal sensitivity (PMS) Q, the fastest apparent convergence (or effective charges) ap- 
proach , nonabelianization methods |^ , Pade approximations |^ , and renormalization-group supplemented Fade 
approximants [||. Of particular relevance is the PMS ansatz which proposes that renormalization-scheme inde- 
^ • pendence, which at lower-loop levels translates into renormalization-scale {^) independence, provides the optimum 
1) \ perturbative prediction of a QCD observable. 

In a particular renormalization scheme (such as MS), the minimal sensitivity principle identifies the appropriate 
choice of renormalization scale jjL for a physical observable as the value at which the observable is independent of 
fjL, providing a method for dealing with the residual renormalization scale dependence that exists in a perturbative 
. calculation truncated to any given order. In explicit calculations, such residual scale dependence decreases as 
' higher-order corrections are included; if a renormalization-group-invariant quantity T is truncated at order (a/Tr)", 
one can easily demonstrate that dF/d/i is order (a/vr)"^ . Consequently, perturbative predictions become virtually 
//-independent at sufficiently high orders. For example, in inclusive semileptonic & — > u decays, a clear progressive 
flattening of the decay rate as a function of /i is observed in the explicit calculations up to two- loop order , and 
this property persists when a Fade estimate of the three-loop correction is included [Q . 

In this paper we demonstrate for semileptonic b ~* u decays that the progressive decrease in renormalization- 
scale dependence with increasing loop-order places bounds on unknown higher-order perturbative coefficients. 
Combined with minimal-sensitivity, the prediction of the decay rate devolving from these bounds leads to an 
estimate of higher-order perturbative corrections to the decay rate. 

The perturbative contributions to the inclusive semileptonic b u decay rate T{b —> uv(t~) may be expressed 
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— r {fi,mb{n),x{fi)) = ml{ji)\^l+ [ao - ai log(w)] x{fi) + [bo - h log(w) + 62log^(w)] x^(^) 
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In the MS scheme in which mf, is identified with a scale-dependent (running) 6-quark mass for four or five active 
flavours, the one- and two-loop order coefficients within (|^) are given by Q 



all Hf : ao — 4.25360 , oi = 5 , 

= 5 : bo^ 26.7848 , bi = 36.9902 , 62 = 17.2917 , 
= 4 : bo = 25.7547 , bi = 38.3935 , 62 = 17.7083 . 

The running mass mb(/i) appearing in (j^) is also known to full four-loop order 

uib ix[fi)\ = nib [x [fio)\ 



c [X i^io)] 

Uf = 4:-. c[x] = a;i2/25 ^ 1.01413a; + 1.38920x2 + 1.09052a;- 
Uf = 5: c[x] = a;i2/23 ^ 1.17549a; + 1.50071x2 0.172486^3 ... ] 
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Renormalization-group invariance of the decay rate F determines a subset of the unknown three-loop Ck coeffi- 
cients and four-loop dk coefficients. From the renormalization-group equation 
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we find the following expressions for terms proportional to x^ and x^ in the final line of 
ci = 26o/3o + ao/3i + 7o (56o - 26i) + 71 (5ao - 2ai) + 572 , 
C2 = i [26i/5o + ai + 571) + 70 (56i - 4&2)] , 

C3 = y (2/3o + 570) , 

di = ^200 + 2/3i5o -I- 3/3oCo -f 70 (5co - 2ci) + 71 (55o - 26i) + 72 (5ao - 2ai) + 573 

1 ^ ^ , 3 ^ „ \ /5, \ 5 

«2 = 7:P2ai -I- Pibi + -Poci -I- 70 ttCi - 2c2 + 71 7701 - 262 + ■;r72ai 



^3 = + -7,^1^2 + 70 77C2 - 2c3 + -7162 



3^ 5 

d4 = |P0C3 + |70C3 



(8) 
(9) 
(10) 
(11) 
(12) 

(13) 
(14) 



2 



Upon substitution of (||) and the four-loop MS results for the /? function ||] and anomalous mass dimension 7 
into (p|-p^, we obtain the following numerical values for these higher- loop coefficients: 



c[^^ = 263.839 



.(5) 
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11.25c, 
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1580.26 , 4"^^ 



765.844 , d' 
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c{' = 249.592 



.(5) _ 



178.755 



.(5) _ 



50.9145 
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10.75c, 
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8.28683 , 



1376.68 , 4^^ 



667.838 , d 
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152.181 



136.833 
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where the superscript denotes the number of active flavours (either = 4 or = 5). In the energy region 
spanning the threshold between = 4 and Uf — b there are four higher-loop parameters that remain undetermined: 

,(4) (_5) ,(4) ,(5)\ 



constant and quark-mass threshold matching conditions imposed at = ||l(](],p] 



dp |. However, continuity of the decay rate at this threshold in conjunction with the coupling- 
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0.633 
-I- 1.9464 
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effectively reduces this set to two unknown parameters. Using the benchmark value (M^) = 0.119 p^ , mb 
4.2 GeV |12[, and the four-loop /3 function, we first find that 



x^''' (4.2 GeV) = 0.07261 , (4.2 GeV) = 0.07269 

ml^^ (4.2 GeV) = 4.208 GeV . 

Imposing continuity of the decay rates to O (x^) at this flavour threshold then yields 
15.66 . 



1.0124^) 



(21) 
(22) 
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Such a continuity condition cannot be extended to O (a;^) , because the first unknown term in the mass threshold 
condition ( pO| ) contributes at O (^x^) to the do coefficient. In the absence of this information, we assume 



d 
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7(5) 
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reflecting the near-equivalence of the known coefficients {d2, da, ^4} for four and five active flavours.^ Thus the RG 
analysis combined with continuity of the decay rate at flavour thresholds implies that the four-loop expression (^ 

effectively has only two undetermined parameters Icq^"*, 4^^| energy range spanning the threshold between 

Uf = A and n/ = 5. 

The decay rate F is a truncated perturbation series which necessarily exhibits residual renormalization scale 
dependence. In general if F is known to 0{x"-), then dF/d/z = 0(0;"+^). Consequently, for any reasonable 
perturbation series, the residual scale dependence diminishes with increasing loop order. A measure of the residual 
scale dependence in the natural energy region^ mb/2 < /i < 2m;, is provided by the deviation of F from its average 
value 



TOfc 



b/2 



(r) 



- 1 



d/i 



(25) 



^In the MS scheme the scale fi = mi, is defined by m;,(mi,) = mj. 

^In (po[), if the coefficient of [x'^' (m^)]* were 20, reflecting a factor often increase between the two previous orders, the difference 

between ^"'^'^ 4^' would be approximately 150, which is small compared to the scales of do we obtain in our analysis below. 

^This energy region is exactly that considered in where the residual scale dependence of the two-loop decay rate was originally 
considered. 
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where the pre-factor of Zm^jj^ leads to a dimensionless x^, and where (F) is the average value of F over the same 
energy interval as (25): 



2rt\u 



(F) = ^ / F(M)dM . (26) 

mi,/2 

The progressive decrease in F's scale dependence implies that must decrease as the loop order of F is increased. 
However, since at three-loop order F (and hence x^) "^iH depend on the parameter c'^\ and at four-loop order will 
depend on |cq'*\ '^o'^^I' the progressive decrease in as loop-order increases necessarily provides constraints on 
these unknown higher-loop parameters. 

To obtain such constraints, we use the central values Ug (Mz) = 0.119 [Q, nib = 4.2 GeV [Q, the four-loop 
(3 and 7 functions, and the threshold matching conditions ([l9[), ( po|) to evaluate a{fi) and nibdi). We then utilize 
discretization to evaluate the integrals and (p^ ) and obtain values shown in Figure 0. The two-loop result 
is by definition independent of c^^ ; the three- loop result depends on Cq^' as displayed in Figure |l|. The requirement 
that the three-loop is less than the two-loop is satisfied provided 

-150 < c^^^ < 290 . (27) 



The Fade estimate Cq*^-* — 188 of ||7[ is well within this interval. Combining the result ( |27D with ( |23[ ) we find the 

(5) 

following constraint on Cq : 

-136 < < 309 , (28) 

(5) 

which, in the notation Cq = 200 -I- A, corresponds to —336 < A < 109 that includes the range of estimates 
— 100 < A < 80 discussed in The optimal value of i.e., the value that minimizes x^i occurs at Cq'''' = 57 
(cg^^ = 74) corresponding most closely to the large-/?o estimate (A = —100) of ||l3|, 

The progressive decrease in x^ to four- loop order restricts Icq^', c^q^'I parameter space to the region indicated 

in Figure ^. Further insight into the allowed d^^^ range can be obtained by employing the Fade estimate c^^ — 188 
01 to determine a constant three-loop value for x^- Upon comparison with the four- loop x^ dependence on the 
unknown parameter dp*' , we find from Figure |^ that x^ continues to decrease with increasing order provided 

-500 < d^a^ < 2500 . (29) 

Although the range of parameter space in Figure § [and in (^7|)-(29)] might superficially be dismissed as too 
large to be of interest, it actually provides valuable control over higher-order corrections in the prediction of the 
actual phenomenological decay rate. If we extract the minimal-sensitivity prediction of the decay rate over the 
entire parameter space of Figure |^, |^ we obtain the following range for the perturbative contributions to the inclusive 
semileptonic decay rate: 

^ ^ 2050 ± 270 GeV^ . (30) 
K 

By contrast, the perturbative two-loop rate, obtained by truncation of the series (0) after its (fully known) O (x^) 
contributions, does not exhibit minimal sensitivity at all in the mb/2 < /i < 2mb region, but is seen to decrease 
with increasing /i from 2053 GeV^ to 1622 GeV^. Indeed, one unanticipated effect of including three- and four- 
loop coeflicients, both known {ci-3, ^2-4} a-nd unknown, is to ensure that minimal sensitivity occurs in the range 
TOb/2 < fi < 2mb over almost all of the parameter space of Figure This property is shown in Figure ^ which 

^Minimal sensitivity occurs when dF/d/j = in the region mi,/2 < < 2mi,. When there exists more than one critical point, we 
choose the one with the smallest second derivative, corresponding to the least sensitive choice. 



4 



compares the two- three- and four-loop contributions to the reduced rate for a typical point in the Fig. |^ parameter 
space. 



The central value of the bound (30) corroborates the Fade estimate 2071 GeV^ obtained in |^. The ±13% 
uncertainty in the rate (|30|), obtained from the combined constraints of a progressive decrease in residual scale 
dependence and minimal sensitivity, is a genuine reflection of the uncertainty following from unknown higher-order 
QCD corrections. Such corrections, which typically are legislated away via series truncation after known-order 
contributions, cannot be disregarded as a source of additional uncertainty in predictions following from perturbative 
quantum field theory. 
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Figure 1: The quantity as a function of Cq (solid curve) for the three-loop contributions to T. The straight 
dashed line represents the two-loop x^, and the criterion that the three- loop is smaller than the two-loop x^ 
constrains Cq^ to the region between the intersection points of the three- and two- loop curves. 
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Figure 2: The dots correspond to points within |cq^\ d'^^^ parameter space for which decreases as loop-order 
increases from two to four. The sharp cutoff in the Cq^' direction corresponds to the range obtained from Figure ^ 
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Figure 3: The quantity as a function of (dashed curve) for the four-loop contributions to F after input of 
the Fade estimate c^q^ = 188 . The straight sohd hne represents the three-loop x^- The criterion that the four- loop 
X^ is smaller three-loop x^ constrains dg^^ to the region between the intersection points of the curves. 
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Figure 4: Two- loop (solid curve), three- loop (dotted curve) and four- loop (dashed-curve) values for the reduced 
rate T/K with c^"^ — 80 (corresponding to a central value A = 114 from |jl^) and d^,*'' — —500 (central value in 
the Fig. 1^ parameter space for the chosen c^^"^ ) . Progressive decrease of residual scale dependence with increasing 
loop order is evident, as is the existence of PMS points for the three- and four-loop curves. Note the absence of a 
PMS point for the two-loop curve. 
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